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Abstract. We show how Z-invariance in the chiral Potts model provides a strategy 
to calculate the pair correlation in the general integrable chiral Potts model using only 
the superintegrable eigenvectors. When the distance between the two spins in the 
correlation function becomes infinite it becomes the square of the order parameter. In 
this way, we show that the spontaneous magnetization can be expressed in terms of 
the inner products of the eigenvectors of the N asymptotically degenerate maximum 
eigenvalues. Using our previous results on these eigenvectors, we are able to obtain the 
order parameter as a sum almost identical to the one given by Baxter. This gives the 
known spontaneous magnetization of the chiral Potts model by an entirely different 
approach. 



1. Introduction 

In 1988 Albertini et al. p!] conjectured a simple formula for the spontaneous 
magnetization M.^. of the chiral Potts model. It took many years to find a proof of this 
conjecture, as the usual corner transfer matrix technique [3] could not be used, because 
the variables on the rapidity lines of the chiral Potts model [21 B] hve on a higher genus 
curve and thus do not satisfy the typical difference property. This conjecture was finally 
proven by Baxter in 2005 j5l [6] using functional equations and the "broken rapidity 
line" technique of Jimbo et al. [7J, invoking two rather mild analyticity assumptions of 
the type commonly used in the field of Yang-Baxter integrable models. Most recently, 
in a series of papers [51 El [101 IHl [121 [13], an algebraic (Ising-like) way of obtaining M.^ 
has been given, providing more insight into the algebraic structure. 

It now looks very plausible that the pair correlation function is also calculable. In 
this paper we shall outline a strategy to attack this problem. As a first application we 
shall derive the order parameter in a new way. 

I Supported in part by the National Science Foundation under grant PHY-07-58139 and by the 
Australian Research Council under Project ID: LX0989627 and DP1096713. 
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1.1. Baxter's Z -invariance for correlation functions 
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Figure 1. Horizontal and vertical rapidity lines and pair correlation functions 
(^A^") = 56"''('?i'92,(73,g4,95,'?6), (sa^c") = Ss"'' (?! : '72, 93, 94, P3, P4, P5 , Pe) and 

(ScSb") = 56™^ (95,96, R.P3,RP4,RP5,RP6 )■ 

The integrable chiral Potts model [3] is Z-invariant by definition — its transfer 
matrices commute with one another and one can permute them without changing the 
partition function Z. Furthermore, Z-invariance also implies that in the thermodynamic 
limit the pair correlation function (s^Sg") in the bulk is given by a set of universal 
functions g^, which can only depend on the 2m rapidities that pass in the same direction 
between the two spins at A and B [HI US] . This may require us to fiip the direction of 
some rapidities by the automorphism q — )■ Rq, given by Xq yq, yq ^ uxq, to make all 
the 2m lines pointing in the same direction. Each spin s can take the N values s = u'^, 
with a = 0, . . . , N — 1, u = e^^^^^ . Shown in figure [l] are three examples with universal 
functions g^'^^ and gg^\ 

The Z-invariance property means that we can calculate the correlation functions of a 
general integrable chiral Potts model on an infinite lattice from special correlations in the 
much simpler superintegrable case [151 [171 HSl [IS] ■ We can take an infinite square lattice 
with special vertical rapidities p and p' alternatingly, but with more general horizontal 
rapidities qi,q2, - ■ ■ , <l2m- Choose two spins within faces in the same vertical column. As 
such a pair correlation function only depends on modulus k and the rapidities passing 
between the two spins, it is independent of p and p', see also the example shown in figure 

gas {s:,s-^-). 

1.2. Superintegrable chiral Potts model 

If we assume alternating vertical rapidities Xpi = yp, ypi = Xp, fip/ = I /ftp, the model 
becomes "superintegrable" [IHl [13, HSl US], as it has Ising-like properties. We shall 
consider here only the case with p = p', so that /Xp = 1. This case both obeys Yang- 
Baxter integrability ^ and has an underlying Onsager algebra PS], whence the name 
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superintegrable [T] . Its spin chain Hamiltonian was first constructed by von Gehlen and 
Rittenberg [18]. 

To calculate pair correlations of the superintegrable spin chain in a ground state, 
we only need correlations in a horizontal row of figure [T| with uniform vertical rapidities 
satisfying Xp = yp. This pair correlation is the expectation value of the spin pair in 
that ground state. From our previous work [201 Ell l22l ES] we know the N ground state 
eigenvectors in the commensurate phase that are needed. 

Several partial results on the pair correlations in the superintegrable quantum spin 
chain already exist, including some results for the leading asymptotic long separation 
behavior from finite-size calculations and conformal field theory f2M 122] , a few terms in 
series expansions in the coupling constant [261 EZl Ell EHl |30] , results using the density 
matrix renormalization group technique [31] and a few exact results for very small chains 
|32j . Using the exact knowledge of the eigenvectors much more can be done for the 
superintegrable spin chain, but in this paper we shall rather go in the vertical direction 
in figure [T| as this will allow us to describe what is needed for the more general integrable 
chiral Potts model and to express the order parameters in terms of inner products of 
ground state vectors. 



1.3. Correlation functions and order parameters 

First we define the row-to-row transfer matrices with periodic boundary conditions in 
the usual way [1], i.e. 

L 

Tg = T{Xg,yg)„^' = jj Wpg{a J " (T'j)W pg{(T J ^-^ - ff'j) , 

J=l 

L 

fr = f{Xr,yr)a'a" = JJ Wpr{(x'j - a"j)Wpr{a'j - (Tj+i), (1.1) 

J=l 

where the chiral Potts Boltzmann weights depend on spin differences modulo N . 

Next, we define the vertical pair correlation functions, between spins Ui and d"i in 
the first column separated by 21 horizontal rapidities, as 

g'^;^[k-q,,---,q^,) = ^Y.^<'^''''^ n Wp,{a-a')Wp,{a-<y"), (1.2) 

{it} all bonds 

where r = l,...,A^ — 1. In the Ising case = 2 we only have r = 1. Assuming periodic 
boundary conditions in both directions, we can write 

^ I M 

9i',\k;qi,---,q2i) = ^Tr ( J] T,,,_,T,,^.) ( H ^.^.-i^...)] , (1-3) 

j=l j=£+l 

M 

with Z = Tr ( I I Tq2^_^Tq.^ \ . Once the thermodynamic limit L, M — )■ oo is taken, 

this should give the most general pair correlations in integrable chiral Potts models on 
infinite Z-invariant lattices. 
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As in [9], the operators Z and X are defined as Z|cr) = oj'^\(y) and X|ct) = and 
L pairs of copies Zj and Xj, acting in the jth column, (j = 1, . . . , L), are introduced. 
They are different from Zj and Xj of [23] which act on the edge variables rij = cTj — aj+i, 
instead of the spins a. The transfer matrices are invariant under the spin shift operator 

L 

X = l[Xj, X\ai,...,aL) = Wi + l,...,aL + l). (1.4) 

i=i 

In fact, the transfer matrix elements only depend on differences Uj = cXj — aj^i in both 
rows and m = ai — (T[ in the first column, i.e. 

( a[, T \a„ {n,}) = Tia,-a[, {n'^}, {n^}) = Tim, {n'^}, {n^}), (1.5) 

writing 

L 

\ai,...,aL) = \{(rj}) = \ai,{nj}), rij = aj-aj+i, ^nj = 0. (1.6) 

i=i 

We define a new (Fourier-transformed) basis [20] 

N-l 

\Q ; {n,}) ^ ^'"^^ 1^1' K»- (1-7) 

0-1=0 

These are eigenvectors of the spin shift operator, as 

ki, {^i}) = ki + l, implies X \Q ; {nj}) = uj'^\Q ; {uj}). (1.8) 

In this new basis, the transfer matrix elements become 



i=Ooi=0 



Af-l Af-1 

or 

N-l 

( Q ; T IQ ; {n^}) = T(m, {n'^}, {n^}) ^ TqHu'^}, {n^}). (1.10) 

m=0 

Thus we find that in the spin-shift Q sector the transfer matrices only depend on bond 
(link, or edge) variables. 
The inverse relation is 

{W-}\T\{a^}) = (1-11) 

Q=0 

Matrix products can also be rewritten in the Q basis: 

N-l 

ctJ'=0 {n"} 
^ N-l N-l N-l 

= ]^EEEE-'""--'''r«(K},K}) 

a'^=0 {n"} Q=0 Q'=0 

X a;Q'K-")fQ,(K},{n^.}) 
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7V-1 N-1 



At E E E ^'"'-"^-"'^^QMin'^}. K})%(K}, {n,}) 

{nil} Q=0 Q'=0 

1 



{„"} Q=0 



l5^u;«K-i)({n;.}|rQfQ|{n^.}). 



:i.l2) 



Q=0 



The pair correlation function can be worked out similarly. In the special case of equal 
horizontal rapidities (g^ = q) and using ({crjlZfllcr}) = w'''^^, we find 

1 _ ^ Af 



92e ik;q,--- ,q) = ^ Tr 

Z {cr} 



j=i i=£+i 



92j 



-V Tr 



Af- 



:i.i3) 



with P = Q — r mod A^. Let the eigenvectors of the transfer matrices be given by 

TQ{x,,yMx,,y,)\yf) = {Aff\yf), {y?\yf) = (1.14) 

where denotes the jth eigenvalue, and let A^^^^^ be the maximum eigenvalue of the 
transfer matrix Tq. Then, in the limit of an infinite number 2M of rows, the partition 
function becomes 

N-l 



^ = E(^-xr'^^«aJ'*^ as L,M^oo, 

Q=0 

as the A^j^^ for <Q < N—1 are asymptotically degenerate as L — )■ oo. Therefore 



1.15) 



N-l J 



92eik;q,---,q) = ^EE 



=0 j=i 



P -\2l 



AO 

max 



max I ^7 ) I ^ max) ) 



:i.i6) 



where J = \ L oo. 

In the limit £ — )■ oo, the pair correlation becomes the product of order parameters. 



N~l 



"^{^-'n = ^ E(^-x|3^max)(3^max|3^Sax), P ^ Q ' T mod N . 



[1.17] 



For the Ising case (A^ = 2) this reduces to the simple formula for the spontaneous 
magnetization rugp = |(3^maxl3^max)l the ordered phase. 



1.4- Baxter's approach 

Recognizing that the A^-state superintegrable model resembles the 2d Ising model 
[U [TTl EH [23l |33l IMl EH] , and that the two parts of its Hamiltonian generate the Onsager 
algebra [HI |33] , Baxter concluded that there must an algebraic route to calculate the 
order parameter Air- He started [8] with a new algebraic approach to calculate J^r 
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for the Ising model, representing it first as a square root of an L by L determinant, 
whicli lie then reduced to an m by m determinant [8], where m < L/2. Having these 
new Ising results Baxter tried to guess their extensions to the iV-state superintegrable 
model. Thus he conjectured [9l[T0] that the special combination Dpq = AiriZpEgY^"^, 
with Zq denoting the partition function with fixed boundary condition of type Q, is 
both a determinant and also given by the sum 



I s\ I sU I s'l 

yi z/2 ■■■ym' ■ 



:i.i8) 



Here m = nip and m' = mq, with mq = [{N — 1)L/N — Q/N\ and mp given similarly. 



and 



S — {Si, S2, . . . , Sra}-, 



0,1, 



:i.i9) 



4) 



= n -^;) = n ik^^ 
= niKs-^.) = n n 

i&w j&v i ev j aw 

in which for a given set s, V" denotes the set of integers i such that Si 
such that Si = 1, while V, W are defined similarly for the set s'. In (1.18) the number 



[1.20) 



and W the set 



of elements in W equals the number in W . The variables q are related to the roots Wi 
of Pp{w) and c[ to the roots w[ of Pq{w), with both polynomials defined by 

iV-l r /. .N\ n L 



E 

n=0 



t 



tun 



:i.2i) 



and the and c'^ by the relation c = {w + l)/{w — 1). This form (1.18) was recently 
proven by lorgov et al. [13j. Finally, in [12], Baxter showed that this sum Dpq is the 
conjectured determinant, which he had already evaluated in the thermodynamic limit 
in [Tl]. Thus he completed an algebraic proof of the order parameter. 



1.5. Approach of lorgov et al. 

Very recently lorgov et al. [13] derived matrix elements of the spin operator in the finite 
superintegrable chiral Potts quantum chain in factorized form. They also gave another 
method to derive the order parameters in the thermodynamic limit, performing the sum 
without going to the determinant formula of Baxter. 



2. Our Approach 

In our previous papers [2n 123] we have given the 2™'3 eigenvector pairs in the ground 
state sectors of the transfer matrices Tq and Tq. Using the notation of Baxter [TUl [TT] . 
but without going to Baxter's reduced representation of the vector space, we label these 



Order Parameter in Chiral Potts Model 7 

eigenvectors and with s' chosen as in (1.19). In fact, the eigenvectors are 

given in [211 123] as 



mq 



j=i i&w;, j=i iGWi, 

where is a subset of {1, 2, ■ ■ ■ , mg} containing n integers. They satisfy the eigenvalue 
equations 

rQ{x,,y,)\X?) = Aj|3^J), fQ{y,,x,)\y^) = A^,\X?). (2.2) 
Particularly, we shall use the ground state eigenvectors 

mp mp 

for which W(d = is the empty set. The "rotation" S = Ylj'^j given by (IV. 124) or 
(IV.150)|with 

S,,P = ^(^if + 4f )1 + K^if - )H,,P + .if Ejp + 4f ET^ , (2.4) 

.if )H,, p - sif E+p - sif E-p , (2.5) 



j,P ~ 2V''22 "T .11 



where E^ p and iij,p are generators of SI2 algebras. Furthermore, we have shown in 
[SI ES] that for the state \Q) = \{nj = 0}) we have E" gl^) = 0, H^,q|1^) = 
(n|E+Q = 0, (^|H,,Q = -(^]| and 



{n\E 



£,0 



a: 



J2 ({^.}k^^^"^GQ({n,},^,,Q) 



{0<nj<iV-l} 
T^lH h'ti'ij—N 



A 



J2 ^-^^'^'GQi{n,},z,^Q)\{n,}), 



(2.6) 



{0<7lj<iV-l} 

nj^H ^nj^—N 



see (IV.66) and (IV. 67). The polynomials here are given in (III. 16) as 



m — 1 



m—1 



GQ{{n,}, z) = J2 KnN+Q{{n,})z'', Gp{{n,}, z) =J2 i^n^+p(K})^", (2.7) 



n=0 



71=0 



with m = mp, m' = mq and coefficients given by (III. 7) and (III. 8) as the sums 



K^{{n,}) = E n 

{0<n^.<JV-l} j = l 
n'j^H yri'j^=l 

u{n,}) = E n 

{0<n^.<JV-l} j = l 
n'j^H hn^=« 



n'- 



i-i 
1=1 



l=j+l 



(2. 



II Equations in our papers [2Ql [21] [221 ES] are denoted here by respectively prefacing I, II, III, or IV to 
the equation number. 
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2. 1 . Form factors 

From (Q, and (|2^ we obtain 



mp 



m' m 



= n ^ n n 11(1 + ^etq) n(i + ^.ejp)!^^), 

where m = mp, m' = nig, and 



u 



^12 _ ^12 _ -^12 „./ _ ^1 _ '='21 ^"^ 



^ _12 

^22 ^22 ^22 



21 



11 ■^ll 11 



(2.9) 



(2.10) 



see (II.C.2), (II.C.IO) and (IV.152). 



Next we expand the products in (2.9). Since 

k 

^E+p|^])Ge{|K•}), Y.^, = kN] 

s=l j 

I 



s=l j 

we find that the only non- vanishing terms in this expansion are those with equal numbers 
of creation and annihilation operators. Consequently, 



m m 



j=i t=i 



l<ii<--<i„<m l<^i<---<^„<m' 



s=l s=l 
m! m 



+ --- + {u,---uj{u\-.. u'^m n ETq W E+p|fi) 

j=i ^=1 



with 



m m 



i=i i=i 



'22 



(2.12) 



(2.13) 



The last term in (2.12) is identically zero unless m = m' . Let Ap = /i^ = 1 in (III. 149) 
and (II. C. 5), and denote = Zjp, Zj = Zj^q, 9j = 6^p and 6'^ = 9jq as in Baxter's 
papers. Then (2.10) becomes 

1-k' 



z^il-k') 



(2.14) 



compare (III. 124) and (III. 126) with Ap = 1. The 6j in these equations are not the same 
as those in Baxter's papers OHO], but are related to our Zj by (II. C. 5), i.e. 

^29, ^ ^-29, ^ ^/ ^ ^/^/ _ _ k'fzj/k'. (2.15) 
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Solving for e^^^ while using Baxter's notations (BaxII.3.16) and (BaxII.2.18)[^ 



1 + A? - (1 + k 
X = l + k'^ + 2k'-—^ or zj = -^ 



/\2 



'1 - k 



i\2 ■ 



(2.16) 



identifying Wj = 1/zj and Cj = cos(^^^) 
9j in Baxter's paper), we have 



k' 



A,-l 



k' 



'{1 + Zj) / {1 — Zj), (where we use 6^ to denote 



(2.17) 



Consequently (2.14) becomes 

l + k' 



A' - l-k' 



A^. + 1 + fc' ' 



1 - k' 



(2.18) 



where the second equation in (2.16) is used to get rid of the factor z'-. Comparing with 
(BaxIV.4.13) and (BaxIV.4.14), we find u'^ — )■ —yj and Uj — )■ —y'y The minus signs 
cancel out upon multiplication. We shall show later that this gives the correct result, 
due to the identification Zj = l/wj. 



Adopting the notations of Baxter as shown in (1.19) and (1.20), we find Wn 
'1? ■ ■ ■ y^n} and = {ji, ■ ■ ■ ,jn}, so that (2.12) becomes 



■■■u: 



X 



n n Hpi^wm'^ 

Similarly, for P -H- Q, we have 



■ ■ - u: 



(2.19) 



X 



(^1 n n Kq\^)^i''^'2'' 



where, instead of (2.10) and (2.13), 



JO 
'22- 



(2.20) 

(2.21) 
(2.22) 



It is easily seen from ( [2lo| ), ( |2l4l ) and ( [2^ , followed by ( [2l8| and ( |2l6| , that 



A. 



1 



A,. + l 



k' 

~k'' 



Finally from (11. C. 6) and (2.17), we have 

(A, + 1)^ 



3 3 
■511-522 



g26»j 



k' 



k'^ 



U 



a;. - 1 + k' 



2 sinh 26 j 



AX, 



(2.23) 



(2.24) 



% Equations in Baxter's papers [SI [9l [101 [H] [El are denoted here by respectively prefacing BaxI, BaxII, 
BaxIII, BaxIV or BaxV to the equation number. Our transfer matrices and X are transposes of the ones 
of Baxter, see footnote to (IV. 107). Therefore, to compare with Baxter, we must replace Q ^ N — Q 
or Zj ^ 1/ Zj = Wj . 
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which relates to the inverse of Zp in (BaxIV.3.7), as it should. 



Looking at (2.19) we now need to evaluate the expectation value in the state \Vl) 



of the product of n creation operators E^p and n annihilation operators E^ g, which all 
operate on the edge variables. 



3. Proposition 



Comparing (2.19) with the sum Dpq of Baxter [HI HO] given here in (1.18), we propose 

(3.1) 



the following identity 



Cf:Df.l 



where A, B, C and D are defined in (1.20) [or (Baxlll.3.44)], while 

As,s'= n 11(^ -4) ^s,s'= n n 



[z'i-z'j). 



(3.2) 



B, = nn 

The difference between A and D of Baxter and A and D here is the replacement of 
c = {z + l)/{z — 1) by z. In B and C, we have also flipped the signs to make them more 
symmetric. 



We shall first consider the simplest case, with n = 1, then prove (3.1) by induction 



3.1. Proof of (3.1) for n = 1 



From (2.6), we find 



(o|ETQE+p|^]) 



^*-0^*-0 {0<nj<JV-l} 
njH t-n^^iV 



Yl GQiM^^j,Q)Gpi{ni} 



Similar to (11.72) or (III. 45), we introduce the polynomial 

h«'^(4, z)^ Yl GqHu,}, zi)Gp{{n,}, z). 



(3.3) 



(3.4) 



{0<ni<JV-l} 
njH hni=JV 



Substituting (2.7) into this equation, we find 

m' — l m—1 

where 



E Ke{{n.})K,{{n,}). 



(3.5) 



(3.6) 



{0<rn<N-l} 
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We can generalize the identity (III. 37) or (III. 44) to P ^ Q cases. For P > Q we have 

QeN+Q,jN+P = QjN+P,m+Q 

= E [(^' - ^ + l)A^Af+i+,_„ -in- £)A« (3.7) 

n=0 



which becomes (III. 37) for P = Q after replacing n j —n. The proof of (3.7) is much 
harder, however, and in the Appendix we shall show why. We have been able to prove 
the identity using a different approach, which is presented in detail elsewhere [36] . Due 



to the symmetry Qi^k = Qk/ given in (3.7), we find 

m' — l m— 1 

^^'"'(4, ^) = E E 4' ^'S,^^P,m+Q = 4)- (3.8) 



£=0 j=0 



Inserting (3.7), and interchanging the order of summations over j and n, we rewrite 
(3.5) as 

m' — l m—l m—1 

^^■"^(4, 4 = E E E 4'^^[(J -n+ l)A?Af+i+,-„ -{n- i)A?^,^,^r.K] 

e=0 n=0 j=n 

m'—l I I m—l m+l—n 

= E (7) E E [(^ - - (- - ^)aM.'^+-\ (3.9) 

1=0 n=0 i=l+l 

where the summation over j is changed ioi = j — n + i+1. Since m' = m or m' = m + 1, 
and also = for n > m' and A^ = for > m, we may extend the intervals of 
summation to < £, n < m', (also noting Yl\=k '^i = ~ Yl'i=l+i I < k), so that 



H«'^(4,4 = EE E (7) -^^"-'[(^-^)An^Af- (n-^)AM. (3.10) 

e=0 n=0 i=l+l 

To evaluate this we enlarge the summation interval for i from £ + 1 <i<m + ^ — nio 
< i < m' and subtract the contributions of < i < ^ and m + £ — n + 1 < i < m'. 
More precisely, defining 



E E [(^ - ^)A^Af - (n - i)AfA^]z^^--\ (3.11) 

n=0 i=0 

m' e 

/5(4 = E E [(^ - - (^ - ^)A?A^]^^+"-\ (3.12) 



n=0 i=0 



and 



^(4 = E E [(^ - ^)A?Af - (n - £)Af A^K+"-\ (3.13) 



n=0 i=m+£— n+1 



we have 



h«-^(4,z) = E (7) [«(4 -/3(4 -7(4]. (3.14) 

i=0 
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It is easily seen that a{z) = by interchanging n f-)- i in the sum of the second term 
of the summand. Similarly, as in [22], we can show that for n < £, so that z > m + 1, m', 
the sum over i vanishes. This leaves 

m' m' 

7(^)=E E [(^ - ^)A^Af - (r^ - ^)A?ArK-^"-^ = 0, (3.15) 

n=l+l i=m+i-n+l 

which is identically zero when we interchange the order of the summations over i and 
n for the first term in the summand (while noting again A^, = if m' = m + 1) and 
make the interchange n <H- z for the second term. 
Consequently, the only nonvanishing term is 

m' 



Using 



h«'^(4,.) = -E 



1=0 



n=0 i=0 



-(A?A;-A«Af)$3 



E 



u 



1 UL 

a V ^ 

du ^ 



u 



(3.16) 



(3.17) 



and 



we find 



5^Afn^ = Pp(n), Y.^k'ju^-' = P'^{u\ Pq(4) = 0, 

j=0 i=0 



h«'^(4,.) 



Pp(4). 



Substituting (3.4) and (3.19) into (3.3), we find 
{^^^j,Q^i,p\^n Q p 



From (III. 17) and (IV. 68) we have 



Ao 



k=l 
m' 



A 



n 



a:.!!' 



z — z 



k)i h'jfl 



A' 



n 



fc=i 



Substituting these into (3.20) we obtain 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



mi,KM= n 03( n 



k=l,k^j ^ J 



-2,- - ZV 



k=lkf^e 



K^j ^k) 
[z„ - Z,) 



(3.22) 
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This is exactly the form in (3.1) with W^i = {^} and W[ = {i}. 



Because of the symmetry (3.8), we find from ( 3.3[ ) that 



'LP 



(3.23) 



3.2. Proof hy induction 



Now we prove the Proposition by induction. The idea of the proof is inspired by reading 
paper [12] by Baxter. Denote 

^n(W^n,0 = (^1 n Kp n Hq\^)- (3-24) 
We have so far shown that for = 1, the following holds: 



^n{w^.w:,) = ^sw^,w'^) n ^r' n 4 



'3 ' 



(3.25) 
(3.26) 



where A, 5, C and D are given in ( |3.2[ ). Substituting these we can summarize (3.24) as 

n n - 4) n n - n ^i^^^^) 



choosing (/)(z) = 1 for ^^^^^ ^L) and (t){z) = 1/z for ^„(W„, PF^- 

Assuming (3.1), or more explicitly (3.27), holds up to n, we shall prove that it also 
holds for n + 1. 



3.3. Proof for n + 1 
Consider 





n Q 


F." F+ 

^e,Q^k,P 


n 




(3.28) 















with = k} and ly^+i = {ly^, i}, so that K+i = Vj{k} and K+i = K/{n, 

where V/U = {j E V, j ^ U}. The operators E^g and E~g were defined for Q = in 
pi] and for general Q in [23]. For fixed Q they satisfy the commutation relations in 
(11.15) and generate a direct sum of spin 1/2 representations of 5l2 algebras, see also the 
appendix of [2T] . 

Because this implies (E^p)^ = and (E~g)^ = 0, we find that 

V^„+i(iy„+i,C+i)=0 if or keW^. (3.29) 
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From (2.6), we see that to each of the roots z[ of the polynomial Pq(z), we associate an 
operator E^q. Thus £ G W'^ is equivalent to z'^ = z'j for some j G W^. Similarly, k G Wn 
corresponds to Zk = Zj for some j G Wn- 

If we let EJq — )• E~p, (or Q = P, i = i), then we must have z'^ = z^. There are three 
possibilities: If i ^ Wn+i, or equivalently i G Ki+i, then ipn+i(}Vn+i,W^^-^^) = 0, due to 

5f^ jii^ p. If, however E^^ E^p making z'g = z,^ 

|fi)=^„(iy,,0, (3.30) 



the commutation relation [E^p, E^- p] 
then 





Efe,pE^P 











as H^p|fi) = — The third case i = j E Wn is similar making z'^^ = Zj leading to a 
reduction to ^JW^^,!^;:), with = W^_,J{e}. 

Likewise, if E+p ^ E+q, or z,^ = z[, for i G Kl+i, then V^„+i(iy,+i, I^^+i) = 0. 
For Zj^ = z\ we have Vn+i(W^n+i> ^^n+i) = V'nl^^n, whereas for i = j G VT^ we have 
another reduction to i^rSWn-- ^n)? '^i^^ = 

To satisfy all these conditions, we must have 



n 



n 



^1= JEB',. 



n 



n 



n+l 



X 



0(4) 



(3.31) 



where (^i^z) is a function to be determined. Note that (3.31) has all the desired zeroes 
and that all fractions in it equal one for z\ = z^. Because of the reduction to (3.22) for 
n = 1, there can only be simple zeroes. Setting n = 0, noting that '?/'o(05 0) = 1 is the 
normalization of the state we would at first have found the more general 



Zu - z^ 



Zn - Z: 



(3.32) 



instead of (3.31) for this case. But different reductions from ipn^ with n > 1 down to 
n = lead to consistency conditions that are satisfied when (^{z^^^z'^) = (pi^Zf}) / (j){z'(). 



So far in this subsection, the reasoning works for both ipn and ■?/'„, so that (3.31) is 



valid for both. However, we still have to compare with (3.22) and (3.23). This leads to 
the conclusion that we must have (j){z) = 1 for ipn and (j){z) = l/z for ipn- 

Because (3.27) is assumed to hold for n, we substitute it into the above equation 

{K+i, 



(3.31). Using 



n IK- -4 
n 



{Wj^, k}, we have 

-^^)= n n -4) 
=5)= n n (^.- 



i) n n 



(3.33) 



n + l 
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and we can write six similar relations. Thus we find 



n n (- -4) n n 



n '^(^^•) 



n n n n 



Z,: - Z. 



^/ n '^(^^ 



This then completes the proof by induction, establishing (3.25) and (3.26). 



3.4- Inner products 

Using (|3l|, ( [3^241 ), (|3^ and (|2^, we find ( [2l9| , (|2^ become 



where 



QP — -DpQ, 



So Sm ^S,s'Bs s' ,s[ ,s'o 



CsDs' 



u , 



(3.34) 



(3.35) 
(3.36) 



3.5. Comparison with Baxter's sum 

Since c = {z + l)/{z — 1), we find 

-2(c. - c' 



- 



-z,- — 



-2(c: - 



(c.-i)(c;.-i)' ''^ (c^-i)(c,-i)^ 

and similar expressions for — Zj and z[ — z'^, so that 

(-2)"(™'-")^,,, - (2)"(— 



^.9. .9' - 



and 



n(c,-i)™'-"n(s-ir 



iew j&v 
(-2)"(™'-")Dy 

n(c:-ir'-^n(s-i)' 



(3.37) 



(3.38) 



(3.39) 



For m = m', the Cj — 1 and c'j — 1 factors in (3.39) cancel out. For m' = m + 1, we find 



from ( pA6| ) and ( |2J8| ) that 

2 

Ci — 1 



;i - k'f - \] 



Zj - 1 



2k' 



-2k' 



(a; - 1)^ - k'^ 
2k' ' 



(3.40) 



(3.41) 



Comparing with (BaxIV.A.l), one can see that the results agree. This shows that we 
have obtained the sum in (BaxIII.3.48) by a completely different route. We can now 
finish the calculation of the order parameter following Baxter [HI |12] closely, apart from 
some subtle differences due to different choice of parameters. Therefore, we shall not 
present the details of our calculation in the next subsection. 
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3.6. Cauchy determinant 



In the most recent paper p[2], Baxter has proven that the sum in (1.18) is a determinant 



In the same way, we may write (3.36) as 

det[l„ + YBY'B^l 



D 



(3.42) 

where we may choose to use either the variables Cj, c[, and y'^ of Baxter \10\ or our 



PQ 



variables z^, z'-, and u\. With the latter choice we have matrix elements 



- 



Y 



Y'- 



(3.43) 



with Mj and u[ given in (2.18). The constants /j and in this case are related to the 
Drinfeld polynomials via 



ff 



±1, 



2j Zj ^ 



i=l> jj^i J=l> jj^i 

SO that B is orthogonal in the sense that 



z' - z'- 



Pp{z'^/K 



B^B 



Im/ if m > m', 



BB' 



Im if m < m', 



(3.44) 



(3.45) 



see subsections 6.1 and 6.2 of [H]. The result is not dependent on using c's or z's. 

Again, comparing with Baxter we have to take account of Zj = l/wj, [so that Cj 
still lies in the interval (—1, 1)]. For P = both zj and 1/zj are roots of the polynomial 



Po{z). From (2.6) we can see that to assign 1/zj instead of zj to E^q is merely a change 
of convention. For P > 0, if zj is a root of Pp{z), then l/zj is a root of Pj\f-p{z), and 



vice versa. The particular choice in (2.6) is no longer arbitrary; it is chosen such that 
the E^p are to act on the ground state sector of spin-translation quantum number P. 
Due to this complication, there are some minus sign differences presented below. 



Instead of (BaxIV.4.21), we find for m 

^m,m(A, A') 



m 



D 



PQ 



A. 



n 



(A2,A'2)ll(l + A:' + \)(l-A:' + A9^ 



(3.46) 



where, as in (BaxIV.2.8), 



n 



c, — c„- 



n 



l<i<j<m 



l<i<j<m' 



(3.47) 



nn( 

i=i j=i 



c,- — c,- 



Therefore, using (2.24), (2.23) and (3.35), we find for m = m' 



n 



7^(1 + k') 11 7^(A 



(3.48) 
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where, as in (BaxIV.3.8), we define 

^„)_ffi,(A + A,)/2 



(3.49) 



n7=i (A + a;.)/2 

However, because Zj = 1/wj, instead of (BaxIV.5.8) and (BaxIV.5.12), we liave liere 



7^(l _ k') = (1 - A;')™— 7^(l + k') = (1 + k'Y'^-''^/''. 



(3.50) 
(3.51) 



As a consequence, we fiave from (3.48) the spontaneous magnetization given as 

= (1 _ l^'2Y{N-r)/N'' _ 

For m' = m + 1, we find the same result as (BaxIV.A.17), namely 



D 



A, 



PQ 



■m,m+l 



(A, A') 



A 



m,m+l 



(A2,A'2 



Hit 



+ A,)2 _ fc^ 



/2 ■ 



(3.52) 



(3.53) 



so that 



n{i - k')n{i + k') UT=in\ 



T /2\r(Ar-r)/Af2 



(3.54) 



which is the same as (3.52). Thus from (3.52) and (3.54), we see that the terms in 



the sum given in (1.17) are all equal in the thermodynamic limit. 



4. Summary and Outlook 

In this paper we have discussed two approaches to the pair correlation function in the 



chiral Potts model. In subsection 1.2 we noted that the ground state eigenvectors of |2T 



suffice for the calculation of the pair correlation in the superintegrable chiral Potts chain 
in the commensurate phase. In fact, to calculate its correlation (Z^Z^^j^), we need to 
evaluate 



(^inE7.p(zi 



z,^) n 



(4.55) 



fc=i fc=i 
which we do not know how to evaluate yet for arbitrarily large system size. However, 
for A^ = 3 and L = 3, it yields a result identical to that of Fabricius and McCoy [32j. In 



subsection |1.1| we explained that the knowledge of the transfer matrix eigenvectors of 
the superintegrable chiral Potts should suffice for the calculation of the pair correlation 
in the more general integrable chiral Potts model. 

In our previous papers [211 122] we have determined 2'^'^ eigenvectors for each Q 



sector and they are given here in (2.1). We have learned to calculate inner products of 



such states chosen from different Q sectors by evaluating the inner product of the state 
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with n creation operators E^p acting on \VL) with one with n annihilation operators E^- q 



acting on {VL\ given in (3.1). We have used these inner products to calculate the pair 



correlation function in the large separation limit, which gave us the order parameter. 



This derivation uses (3.7), which in earlier preprint versions of this paper had to be left 
as a conjecture, but finally got proved PB] in a rather lengthy manner. 

However, if we want to go further and study the separation dependence of the 



pair correlation function using (1.16), we will need to know the inner products for all 
eigenvectors {yf) that are not orthogonal to the maximum eigenvectors. One way is to 
first calculate all eigenvectors of the superintegrable T2 model and the related higher- 
spin XXZ-model by Bethe Ansatz working out the prescription of Tarasov [2Zl EH ES] , 
and then to make proper linear combinations in each degenerate eigenvalue space [39] 
as we have done 122] for the ground state sector using Onsager algebra and quantum 
loop algebra. But it will be very hard to work out all necessary details this way. We 
believe that there should be a better way using more of the combinatorial structure of 
the problem, which we may start exploring using the approach proposed in subsection 
1.2| for which we have the needed eigenvectors explicitly. 
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Appendix A. Generating Function Q{t,u) 

The generating functions g and g are defined in (11.62), (111.14) and (111.15) as 

{N~l)L-N 

9{{nj},t)= ^rn(K})r, 
m=0 

{N-1)L-N 

9i{nj},t)= ^U{nA)t"'- (A.l) 

m=0 

It was shown in j22j that they have the simple form 

L 

g{{n,}, t) = r (K), t) = {l- t^)^-i tu'^^y'. (A.2) 

We now define the two-variable generating function 

g{t,u)= Yl 9i{nj},t)g{{nj},u). (A.3) 

{0<nj<N-l} 
ni-\ \-n]^—N 



Order Parameter in Chiral Potts Model 



19 



Substituting (A.l) into the above equation and comparing with (3.6), we find 

Af(r-l) N{r~l) 

g{t,u) = Q^^t'^''^ r = {N- l)L/N. 

1=0 k=0 

From Baxter's paper [TU], we reahzed that the coefficients Qe^k must be related to 



(A.4) 



the coefficients of the Drinfeld polynomials (1.21). However, the method that we used 
in [22] cannot be used, as is shown next. 



Appendix A.l. Difficulties for P ^ Q 



From (3.6) and (2.8|), we may write 

Gm+Q,jN+p = 



EEE n 

{0<\-,n^,n'^<N-l} i=l 













rii 



(A.5) 



where hj = J2i>j \ and N'j = Y.^^. n'-. Let /i. = n ■ + n[, so that ^/ii = (£ + l)A^ + Q. 
Then we find 

L 



UN+Q,jN+P 



EEE n 

{0<A,,ni,Mi<iV-l} i=l 







'rii + Xf 


rii 




Hi 



ni{ai-Ni+bi) 



{0<Mj,Aj<]V-l} 

E Mi=(«+i)iv+g, E A,=jjv+p 



EE In{M;{X.}) 

l<Mj,Aj<]V- 

where = X]i<i /^^ ^'^'^ = ^i<j ^^"^ 

x„({M;{A.})= E n 



(A.6) 



{0<n^<N~l} i = l 
H Vnj^^n 







'rii + Aj" 


rii 







ni(ai-Ni)+nibi 



(A.7) 



see (III.22). Using ql+i = ^ = + 1)A^ + Q and 6o = E Ai = jN + P, we find from 
(III.28), (III.33), (III.34) and (III.35) the relation 

E(-l)"c^^"'X„({/.,};{A.})r 



n=0 



6o 



n=0 



where 



Xn({A.};{M)= E n 



{0<ni<N-l} i = l 
fi-^ _| f-n^ — n 







"rij + jji' 






rii 



ni{bi-Ni)+niai 



(A.9) 



as in (III. 26). By equating the coefficients of on both sides of (A.8), we relate the 



sums In and I„. For P = Q, this process simplifies as {u!2^^t;u)j^ = 1 + in (A.8) 
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For P 7^ Q, it is more complicated, but we still can use (A. 8) to relate Xn with Xm and 
other X„'s with the n < N coefficients calculated from (III. 11). Thus we obtain 



Qm+QjN+p — ~ j)A2(.iA^ + Qf^ 



'+1)N+Q,{j-1)N+P ^ ^IJ 1 



where 



E 

A;=0 



N-P + Q 
Q-k 



k^-kP 



(A.IO) 



(A.ll) 



EMj=(<+i)JV+Q+-P-fe 



We have changed the X^s to X„'s by changing the summation variables to Aj = X[ + 
and /Xj = — ra^, where the n/s are the summation variables in (A. 7) and (A. 9). 

For P = Q, the only nonzero term is k = Q = P. Since Xo{{fi^}; {X'j}) = 1, we find 



(A.12) 



{0<m'<JV-1} {0<A'.<]V-1} 



so that (A.IO) becomes identical to (III. 43) as it should be. 

For P > Q, we have to keep using (A. 8) to express each Xfc({/i^}; {A'j}) in (A.ll) in 
terms of Xi{{fii}; {Aj}) with i < k and n = J2i -^j ^ -^i; we have done once before 
to arrive at ( A.10[ ), continuing until n < k, such that Xk = 0. Let 

J2 ^=Cm, CnN+Q = A'^. (A.13) 



{0<Ai<iV-l} 



Since Xq = 1, we find it is possible to express U^j^ as a sum of C(^+i)7v+Q+p-mCjAf-m- 
However, this is very tedious and messy. Using Maple we have discovered (and checked 
for different P, Q and values and small system sizes L) that the following equation 
holds: 



Q,P 



(A. 14) 



n=0 



n=0 



SO that 

i j 
Gm+Q,N+P = £(^ - J + 2n)Ag.,^„Af_„ + J^Wf+r,. 



n=0 



n=0 



j-n-1 



m 



n=0 m,=0 m=0 

E (£ + J - 2m) A,«+i^^._„A,^ + (j - m + 1)A« Af+i+^._^ - (j - m) A,«+,^^_„A 



m=0 



m=0 



5^ (£ - m)A?^,^^._A^ + (j - m + l)A^Af^, 



+j—m 



(A.15) 



which is the identity in (3.7). In the first preprint version |10] of this work (A. 14) had 
to be left as a conjecture. In the meantime, a lengthy proof of (A.15) has been found, 
which is presented elsewhere 
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